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, Using exact diagonalisation and Density Matrix Renormalisation group (DMRG) approach we 

analyse the transition to a localised state of a weakly interacting quasi- ID Bose gas subjected to 
^ ' a quasiperiodic potential. The analysis is performed by calculating the superfiuid fraction, density 

profile, momentum distribution and visibility for different periodicities of the second lattice and in 
the presence (or not) of a weak repulsive interaction. It is shown that the transition is sharper 
towards the maximally incommensurate ratio between the two lattice periodicities, and shifted to 
higher values of the second lattice strength by weak repulsive interactions. We also relate our results 
to recent experiments. 
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C/2 ' I. INTRODUCTION 

CD ' 

Ultracold atoms in optical lattices [l[ represent a powerful and robust tool for simulating simple quantum systems 
with a broad tunability of the parameters, thus serving as quantum simulators[2j to reproduce the fundamental 
phenomena of condensed matter systems. One of the most fascinating and still long-standing problem in condensed 
matter physics is represented by the localisation of waves in disordered media;3]. Actually, evidence of the Anderson 
localisation for light waves in disordered media has been provided by the observation of the modification of the diffusive 
regime, indicating a transition from a conductor to an insulator transitional^]. However, a careful analysis is necessary 
in order to separate localisation from absorption effects [6|. Anderson localisation in three-dimensions has also been 
recently observed with ultrasound waves [7|. In the case of electronic waves in solids, the situation is made more 
complicated by the long range coulombic repulsion between electrons^. Localisation problems can also be considered 
I | in the case of bosonic systems in a disordered medium such as 4 He in porous material^ or disordered Josephson 
■ Junction arrayspLQj. In both of these systems, the bosons are again strongly interacting. Thus, the interplay between 
^\ ' disorder and interaction is an important issue in theoretical condensed matter physics. The combination of ultracold 
r**"" . atoms and optical potentials is offering a novel platform to study disorder and interaction related phenomena, in 
which most of the relevant physical parameters, including interaction, are t unable [111] . Moreover, a disordered optical 
potential can be realized either with laser speckles[l2j or with quasi-periodic optical lattices. [lj]. For both types of 
' disorder, Anderson localisation has been recently experimentally observed, [lil Il5j 

In this work we focus on quasiperiodic potentials, which exhibit properties common to both periodic and disordered 
systems jl6|. As in the case of periodic lattices, their spectrum shows reminiscences of energy bands. On the other 
hand, owing to the lack of any translational invariance, they support the existence of localised states, which behave very 



similarly to the ones generated by truly random potentials 17] . Localisation of noninteracting particles in q uasiperiodic 



potentials has been intensively studied in relation with electronic structures of incommensurate solids. [18j The case of 
a bichromatic lattice of Refs.fll, EH is described by the so called Harper model[l9| or ID tight-binding Aubry- Andre 
model[2(J Hi], [22j the Hamiltonian of which is: 

_ _ i 

H = -t^2(\n)(n+l\ + \n+l)(n\) + A^2cos(2nnf3 + <j))\n)(n\. (1) 

n n 

Here t is the transfer integral between nearest neighbour sites of the main lattice, A is the amplitude of the quasiperi- 
odic modulation of the potential energy, while f3 is an irrational number and <fi is a phase. The Aubry- Andre model 
can be experimentally realized by superimposing two optical lattices of incommensurate wavelengths, the primary 
lattice height, s\, being much larger than the secondary lattice height s 2 - 

V(x) = StEm cos 2 (kxx) + s 2 E R2 cos 2 (k 2 x) (2) 

where k\,2 is the lattice wavevector and E Ri = h 2 k\/2m the respective recoil energy. At variance with the completely 
random potential which localises all the states as soon as it is switched on, [23| in the Aubry- Andre case the strength 
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of the potential must exceed the critical value A/t = 2 to induce a localised ground state with a localisation length 
l/ln(A/2i) in units of lattice spacing. Moreover, the localised state at A/t > 2 is related to the dclocaliscd state at 
A/t < 2 by a duality transformation. The transition point is self-dual, and its spectrum is singular continuous, with 
wavefunctions having a fractal-like structure [20j. 

In the recent experiment in Florence by Roati et al. [TH , Anderson localisation has been observed for a noninteracting 
39 K BEC in an incommensurate bichromatic lattice where the effect of interactions have been cancelled by tuning a 
static magnetic field in proximity of a Feshbach resonance to set the scattering length to zero. The crossover between 
extended to localised states has been studied in detail by looking at the expansion of the BEC and by studying spatial 
and momentum distribution of the atoms, with results in qualitative agreement with the Aubry- Andre predictions. 
In particular, they have verified the scaling behaviour of the critical disorder strength and shown that while the 
transition to localisation is expected to be sharp in the case of a maximally incommensurate ratio (5 = (\/5— l)/2, it 
is broadened for the experimental parameter (3 = 1.1972. 

So far, the incommensurate lattice supports a localisation transition for noninteracting particles, but a richer physics 
arises when interactions are present. A first interesting question is how the critical disorder strength would be modified 
in the presence of a weak interaction. Another question is the nature of the localised phase of interacting bosons. 
Indeed, the repulsion between bosons disfavours the condensation of all particles in the lowest localised eigenstate, and 
should instead favour a phases in which although the bosons are localised, their density is homogeneous. This problem 
has been discussed in the case of a weakly interacting system in Ref. [2J] , where the ground state of the system has 
been calculated for a three-colour optical lattice. Introducing repulsive interactions between the atoms, the numerical 
integration of the ID time dependent Gross-Pitaevskii equation shows that the ground state wavefunction becomes a 
superposition of many single-particle localised states, which add up to form an overall extended state. Similar results 
for a bichromatic lattice have been presented in the work by Lye et al.fl3|. and Quantum-Monte Carlo calculations 
for the case of large interaction strengths have been performed by Roscilde [Ifl [26[ . 

Besides the modification of the critical disorder strength, the combination of lattice and interaction can lead to new 
phases that compete with the localised state. In absence of disorder, bosons on a lattice with repulsive interactions 
display, for commensurate filling, a superfluid (SF) to Mott insulator (MI) transition as the repulsion is increased 271 ] . 
with the superfluid phase displaying large density fluctuations and a gapless excitation spectrum, while the Mott 
phase is incompressible and has a gap in the excitation spectrum. If one considers both repulsive interactions and 
disorder, these two effects will compete: while disorder makes the bosons localised, short-range repulsive interaction 
energy increases as the square of boson density and hence the total energy of the system is minimised by depleting the 
localised condensate towards a more uniform density distribution. As a result, in a lattice Bose gas with short-range 
interactions a novel Bose-glass (BG) phase, non superfluid yet compressible, emerges between the superfluid and the 
Mott-insulator[27jj. For the specific case of a one-dimensional Bose gas subjected to an uncorrelated disorder (in 
absence of a lattice) , the phase diagram has been obtained by Giamarchi and Schulz [28[ , showing that while for zero 
interactions the system is always localised, for nonzero values of the repulsive interactions a superfluid phase is possible 
at small values of disorder. Ref. [28[ also predicted that the non-superfluid (Bose-glass) phase of an interacting Bose 
gas is expected to differ markedly from the non-interacting Anderson-localised (AG) phase, e.g. the density profile 
of a Bose glass phase is rather uniform, in contrast to the highly inhomogeneous density profile of an ideal Bose 
gas in a disordered potential where all the particles occupy the lowest single-particle localised orbital. The case of 
lattice hard core bosons with nearest neighbour repulsion was considered in [291 ] . It was found that a superfluid dome 
could be present for strong enough nearest neighbour attraction and weak enough disorder. The phase diagram of 
a disordered, interacting lattice Bose gas in one-dimension has been the subject of several numerical investigations 
by exact diagonalisation [3(j, quantum Monte Carlo methods [3ll |32|. strong coupling expansions 33] , and Density- 
Matrix renormalization group approaches (DMRG)f34|, that have established the existence of a Mott insulating phase 
separated from the superfluid phase by a Bose glass phase for disorder not excessively strong. For stronger disorder, 
these numerical studies have established that only the Bose glass and the superfluid is present. Also, the existence of 
a superfluid dome in the phase diagram for a quasiperiodic lattice potential has been obtained in the commensurate 
caselM SS El. 

Here we focus on the quasi-lD Bose gas subjected to a quasiperiodic potential in the presence or absence of weak 
interactions and try to analyse the scaling behaviour of the critical interaction strength by exact diagonalisation or 
DMRG. In particular, we characterise the transition by means of the behaviour of the superfluid density, the density 
profile and the visibility. In the case of open boundary conditions, where the superfluid density is undefined, another 
physical quantity is also introduced for the non-interacting case, i.e. the inverse participation ratio and its results are 
discussed. 
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II. THE MODEL 
A. The Hamiltonian 



We consider a one-dimensional Bose gas at zero temperature subjected to a trichromatic lattice potential as the one 
ofEq.Q 



H= J dxi>l(x) ( -— V 2 + V(x)j Mx) 



+ | / dx$l(x)1>l(x)1n,(x)ih(x), (3) 



where tpb{%) is the bosonic field operator, m is the atomic mass and g = AnTi 2 a s /m is the contact interaction expressed 
in terms of the scattering length a s . When the primary lattice height S\ is large compared to the atomic recoil energy 
Er i and to the height of the secondary lattice S2 , we can map the Hamiltonian on a Bose-Hubbard model [37l |: 

N-x v N 

H = -tJ2( b t b i+i + h.c.) + -J2ni(ni-l) 
»=i i=i 

N N 

- Li)^ni+}^Ajni, (4) 
i=l i=l 

where b\, bi are bosonic field operators on the site i, t is the hopping amplitude, U is the on-site interaction, fx 
is the chemical potential, N is the total number of lattice sites; the parameters U, t are related to those of the 
continuum model ©(e.g. see Refs. [13, l38| ) and their respective dependence upon the recoil energy, lattice depth 
and scattering length can be calculated both analytically and numerically [39(. The effect of the second lattice is to 
induce a modulation of the on-site energies according to Aj = A cos(27r/3z + <f>) , with A oc S2 and <fi being a phase shift 

while in the experiments = k^/ki. By taking the parameters of [40( one can estimate s^/t = 2.5 53. 

In order to characterise the localisation transition, we evaluate the following observables in the case of periodic 

boundary conditions: (i) the superfluid fraction as the response to twisted boundary conditions, f s = ^j-t g g i , 
which on the lattice model reads 

/-5^(*"-*Sfc). m 

where Np is the particle number, Epg C is the ground state energy and E^ p is the ground state energy with boundary 
conditions twisted by an angle [3l|, (ii) the spatial density profile (n?) = (bjbj), and (iii) the momentum distribution, 
given by the Fourier transform of the one-body density matrix 

n{q)=UY,^ q{l ' m)a {b\b m ). (6) 



i ni 

with a — ir/ki being the primary lattice spacing and J\f a normalisation constant. 

B. Numerical methods 

The ground state properties of the interacting Bose gas in the bichromatic lattice have been determined by the 
Density Matrix Renormalization Group (DMRG) method]!^ 511 and exact diagonalisation (ED) method (in the 
noninteracting case). 

Concerning DMRG, it is a quasi-exact numerical technique widely employed for studying strongly correlated systems 
in low dimensions. Based on the renormalization, it finds efficiently the ground state of a relatively large system with 
quite high precision. We have been considering a system with periodic boundary conditions and used first the infinite- 
size algorithm to build the Hamiltonian up to the length L, then we have been resorting to the finite-size algorithm 
to increase the precision within many sweeps. Since the Hilbert space of bosons is infinite, to keep it finite we have 
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been choosing the maximal number n max of boson states approximately of the order 5(n) , where (n) is the average 
number of bosons per site, varying n max between n max — 10 and n max = 25. The number of eigenstates of the 
reduced density matrix are chosen in the range 150 — 250. The number of sites in the DMRG is N — 50. The initial 
phase of the quasi-periodic potential is chosen <j) = tt. 

In the nonintcracting case, U = 0, we have been employing the exact diagonalisation method on a chain of 50 up 
to 1000 sites. The supcrfluid fraction depends on the energy change under twisted-phase boundary conditions. We 
checked that the superfluid fraction computed from ([5]) does not vary strongly at critical points A c /t when small 
angles (say 8 = 0.1 rad) are replaced with large angles (say 8 = tt). Thus, the location of the critical point is not 
sensitive to the choice of the twist angle in Eq. ([5]). This allows us to use anti-periodic boundary conditions (i.e. a 
twist angle 8 = tt) to calculate the superfluid fraction (both in ED and DMRG). We have also analysed the variation 
of the superfluid fraction as a function of the shift 4> of the quasi-periodic potential in the Hamiltonian Q. Under 
some initial phases with fixed length of chain, we found that the critical points are shifted to a value much smaller 
than A/t = 2. This behaviour is associated with the formation of an isolated state below the continuum of extended 
states which is localised near the edges. The origin of this localised state is simply the fact that the on site potential 
A, does not have the periodicity of N, thus causing strong fluctuations of the potential near the edges in the system 
with periodic or aperiodic boundary conditions. In our calculations we choose the initial phase such that no spurious 
localised state is formed near the edges, ensuring that the critical points are close to A c /t = 2. In addition, when we 
perform finite size scaling we consider a sequence of periodic approximants to the potential A^ associated with larger 
and larger ring sizes. In the noninteracting limit we consider also open boundary conditions (OBC). In this case the 
superfluid density ([5]) fails to be a good indicator of the localisation transition, but the inverse participation ratio 41 
can be employed. 



III. ANALYSIS OF THE LOCALISATION TRANSITION 



We proceed here to present the results for the localisation transition occurring at increasing height of the secondary 
lattice, with the purpose of analysing the localisation transition at varying the parameter (3, i.e. the periodicity of 
the secondary lattice. Although mathematically the transition is expected only for irrational values of (3 [12], in the 
numerical calculations (3 is a rational number. In the experiment, (3 is fixed by the ratio of the wavelengths of the two 
laser beams used for creating the optical lattices, which are only known up to a finite number of significant digits, so 
that the rational or irrational character of (3 cannot be decided. In principle, if (3 is rational, the system is in a periodic 
potential, and it thus possesses only extended states forming a finite number of bands for any strength of the potential. 
However, if (3 is given by the ratio of two mutually prime integers, with a large denominator, the periodical behaviour 
can only be observed on a very large length scale. Since, in contrast to the ideal case studied in Mathematics 22], both 
in simulations and in experiments the number of sites of the lattice is finite, the very large lengthscale periodicity 
mentioned above may not be accessible in practice, and the question arises on the occurrence and the nature of the 
Anderson localisation transition (or crossover). More precisely, in the case of A 3> t, if the periodicity of the potential 
is N, degenerate perturbation theory provides an estimate of ~ t(t/ A) N ~ 1 for the bandwidth of the extended states. 
On timescales small with respect to the inverse of the bandwidth, particles will exhibit the appearance of localisation. 
We thus investigate here the Anderson localisation behaviour of a finite lattice and rational f3 parameter, which is 
relevant for the ongoing experiments. 



A. Noninteracting case 



We begin our analysis by considering first the noninteracting limit U — 0, In this case, the ground state of the 
Hamiltonian Eq.Q can be obtained by exact diagonalisation. 



1. The inverse participation ratio 



In the noninteracting regime it is possible to adopt several definition of localisation. One which turns out to be 
particularly convenient for our problem is the inverse participation ratio in the ground state[4l]. This can be expressed 



Po = £K (0) l 4 (7) 
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FIG. 1: (Color online) Inverse participation ratio IPR= l — 8w/N (adimensional) for a noninteracting Bose gas in a quasiperiodic 
potential with (3 = (y5 — l)/2 with open boundary conditions as a function of the number of sites in the chain for various 
values of the height of the second lattice A (in units of t), from bottom the values are A/t=1.9, 1.95, 2, 2.05, and 2.1 . The 
inset shows the same quantity in logarithmic scale. 

where ipf^ is the normalised ground state wavefunction of the noninteracting Hamiltonian (f4]) on the basis |i) of the 
lattice site occupation number. The inverse participation ratio is related to the probability of return to the initial 
point [4l[ and allows to distinguish between extended and localised states: it decays as 1/N where TV is the number 
of lattice sites in the extended regime and goes to a constant value in the localised regime. 

The inverse participation ratio can be extracted from the numerical calculations for a system with open boundary 
conditions by computing the average number fluctuations, 

<^ = ^E[< n ?>-w 2 ]- ( 8 ) 

i 

Indeed, using the expansion of the boson site operators bi on the operators b v corresponding to the eigenvectors ifxf 

of the Hamiltonian, ie hi = J2 V V'l' 7 ^! an d the f ac t that m the ground state all the N particles occupy the state ip^ 
we obtain a relation between the number fluctuations and the inverse participation ratio 

5 W = N(1-J2\4 0) n (9) 

i 

Since the inverse participation ratio has two distinct large- N behaviour in the extended and localised regimes, it is 
a good indicator to identify the transition from the extended regime, where Sw/N = 1 + 0(1 /N), to the localised 
regime where Sw/N = 1 — const. FigfT] shows the change in behaviour of the inverse participation ratio across the 
localisation threshold. 

Figure [2] shows the average number fluctuation as a function of the height of the secondary lattice A for different 
choices of the quasiperiodic potential. In the first panel, we choose for /3 subsequent approximants of the irrational 
number (1 + Vo)/2, according to the Fibonacci sequence: Ft = 1, F% = 1, F n+1 — F n + F„_i, as lim^oo F n+ i/F n — 
(1 + \/b)/2. The figure shows that the transition at A/t — 2 becomes more and more pronounced as the order of 
the approximation increases. In the second panel, we show the behaviour of the inverse participation ratio for three 
choices of /? which we will consider throughout the paper: the maximally irrational value (3 = (v5 — l)/2, the value 
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FIG. 2: (Color online) Average number fluctuation Sw (adimensional) for a noninteracting Bose gas in a quasiperiodic potential 
with open boundary conditions as a function of the height of the second lattice A (in units of i) for a number of lattice sites 
N = 200, at varying choices of f3. First panel, subsequent approximants of the number (3 = (1 + VE)/2 according to the 
Fibonacci sequence. Second panel p = (V§ - l)/2 (blue solid line), (3 = 830/1076 ( green dotted line), and (3 = 1032/862 (red 
dash-dotted line). 

adopted in the first Florence experiment [l3| P = 830/1076, and the value used in the second Florence experiment 
PJ] (3 = 1032/862. 

2. The superfluid fraction 

Since the inverse participation ratio is defined only in the noninteracting limit, we further characterise the transition 
by monitoring the superfluid fraction, which is nonzero in the extended phase and vanishes in the localised phase. We 
obtain the superfluid fraction from the change of the ground state energy of the system to twisted boundary conditions, 
Eq.([5]). Figure [3] shows the behaviour of the superfluid fraction as a function of the height of the secondary lattice A, 
for various choices of the quasiperiodicity parameter (3 and for a number of lattice sites equal to 50, which is close to 
the experimental conditions. We have considered three possible values of (3: the irrational number (3 = (\/5 — l)/2, 
and the values (3 — 830/1076, j3 = 1032/862 adopted in the Florence experiments [l3|,[l5[. From the data we see that 
the superfluid fraction vanishes around the value A/t = 2, as predicted by the Aubry- Andre model [2(J, the details 
however depend on the choice of (3. We remark that the transition is particularly smooth in the case of the second 
Florence experiment, where the potential has five minima very close in energy. 

The Anderson localisation transition can be also seen by looking at the spatial density profiles, illustrated in Fig. 
2] at increasing values of A. In correspondence to the transition the profile displays a clear change of behaviour as all 
the particles tend to cluster around the localised state. The position of the Anderson-localised peak corresponds to 
the position of the absolute minimum of the quasiperiodic potential. 

3. Finite size scaling 

All the results of the previous section were obtained for a system of finite length. In such a system, there is no 
actual localisation transition as the localisation length can never exceed the length of the system and the spectrum 
always remains a pure point one. In order to interpolate the results to the case of an infinite system which possesses 
a true localisation transition between localised and extended state, we make a finite-size scaling Ansatz.[44{ We thus 
assume that the superfluid fraction is of the form: 

f s (N,A)=N- a <S>(-^L.y (10) 

where £(A) is a characteristic length that diverges at the localisation transition as |A — Ad" 1 '. The function $(a;) is 
assumed to be regular for i«l and to behave as Q(x) ~ x a for x 3> 1 on the extended side, so that the superfluid 
fraction becomes independent of system size in the thermodynamic limit. 
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FIG. 3: Superfluid fraction for a noninteracting Bose gas in a quasiperiodic potential as a function of the height of the second 
lattice A in units of the hopping strength t for three choices of the quasiperiodicity period j3: maximally irrational value 
(3 — (\/5 — l)/2 (solid lines), first Florence experiment [l3T ] j3 = 830/1076 (dashed lines), and second Florence experiment [l5j] 
(3 — 1032/862 (dot-dashed lines). All the curves are drawn for cj> — it and N = 50. 





FIG. 4: (Colour online) Contour plot of the density profile for a noninteracting Bose gas in a quasiperiodic potential with 
periodic boundary conditions at increasing height of the second lattice A for a number of lattice sites N = 50, at varying 
choices of p: (3 = (VE - l)/2 (a), (3 = 830/1076 (b), and f3 = 1032/862 (c).In the three cases localisation is clearly visible for 
A > 2t. The value of the phase shift has been chosen <f> = 37r/10 for (a) and (b), and <j> = —2ir/5 for (c). 



When approaching the transition from the localised side, this characteristic length is obviously the localisation 
length. On the delocalised state, the identification of the characteristic length is less straightforward. Using the 
Aubry- Andre duality, there is now a dimcnsionlcss localisation length £ defined in momentum space, which measures 
the number of plane waves of quasi-momentum k + n(3 that must be combined in the repeated zone scheme to form 
an extended wavefunction. When these quasi-momenta are folded back in the first Brillouin zone, we end up with a 
wavefunction formed by linear combination of I distinct plane waves, whose momenta are in the interval [—w/a, it /a]. 
The typical distance between these different momenta can be estimated to be of the order of 2ir/ (la). Thus, in order to 
resolve the difference between these momenta, it is necessary to be in a system of size N such that 2n/ (Not) <C 27r/ (lot) 
i. e. N ^S> I. Thus, the length I has the interpretation of the minimal system length necessary to distinguish a extended 
state from a localised state. We thus identify the length lot to £(A)a in the extended regime. Using the exact result 
on the localisation length in the Auhry-Andre model, we have £(A) = 1/| ln(A/2f)| ~ 1/|A — 2t|. Thus, the exponent 
v = 1 on both sides of the transition. Using this expression for £(A), we have been able to collapse our numerical data 
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FIG. 5: Double logarithmic scaling plot of the superfluid fraction in the case of non-interacting bosons (/3 = (y5 — l)/2). The 
data for different system sizes are collapsing on a single curve once properly scaled. Deviations from scaling far away from the 
transition point are visible on the graph, indicating that the corresponding points are already outside of the scaling region. 

on a single curve, as can be seen on Fig. [5] using an exponent a — 1/2. This leads to a superfluid fraction vanishing 
as \2t — A I 1 / 2 as the transition is approached from the extended side in the thermodynamic limit. 

B. Weakly interacting case 

Using DMRG we can access the regime of arbitrary interaction strength (36J. In particular, we consider here the 
effect of weak repulsive interactions. In this regime, interaction reduces localisation, as repulsion keeps particles from 
coming on top of each other. Figure [5] shows the superfluid fraction as a function of the height of the secondary lattice 
A. The comparison with the noninteracting case shows that in the interacting case the localisation transition is shifted 
to higher values of the secondary lattice A. The delocalised phase, having both non- vanishing compressibility and 
non- vanishing superfluid fraction is in a Luttinger liquid phase. [36l l4f| In particular, it displays a quasi- long-range 
superfluid order, with the single particle density matrix decaying as a power law as a function of distance. This also 
gives a power law divergence of the momentum distribution at low momentum in the infinite-system limit [36]. The 
observation of the shift of the transition could provide a sensitive measure of the interaction strength U , hence of the 
scattering length a s . 



1. Finite size scaling 

In the interacting case, we can also develop a finite size scaling Ansatz analogous to (fTTj|) . We now assume that 
£(A) ~ |A — A C (U)\~ U , where A c ({7) is the transition point in the presence of interaction U. The main difficulty 
in applying finite size scaling to the transition in the interacting case is that in contrast to the non-interacting case 
of section UlI A 31 we do not have a priori knowledge of A C (U). However, we note that in the localized phase, f s is 
decreasing exponentially with the system size, and at the critical point, according to finite-size scaling is decreasing 
as a power law. By plotting f s as a function of A with U = 0.02i, we find that f s is a non-decreasing function of N 
for A/t < 3.2. Thus, the localized phase is obtained for A > 3.2t. We have assumed that A c = 3.2t was the location 
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FIG. 6: Comparison of the superfluid fraction for a noninteracting (straight lines, U — 0) and weakly interacting Bose gas 
(lines with circles, U/t — 0.02) in a quasiperiodic potential as a function of the height of the secondary lattice A/t at varying 
choices of the periodicity parameter: solid lines, (3 — (v5 — l)/2; dot-dashed lines, (3 = 1032/862. Note the enhancement of 
the superfluid fraction by the interaction for both cases. For all the curves we have used a number of lattice sites N = 50 and 
a filling factor 1/2. 

of the critical point, and fitted the N dependence of / S (A C , N) to obtain the exponent a = 0.38. We have been able 
to collapse the curves (see Fig. [7]) using v = 1.25. The collapse of the data points is a verification of the scaling 
hypothesis. 

IV. MOMENTUM DISTRIBUTION 

The momentum distribution is one of the experimentally accessible observables [l5l ]. and allows to distinguish 
among the extended and localised regimes. We have calculated the momentum distribution for various choices of the 
quasiperiodicity parameter [3, both in the noninteracting and in the weakly interacting case. 

Our results are reported in Fig. [8] In the extended regime, in addition to the central peak at q = side peaks in 
the momentum distribution appear, at a position related to the value of the periodicity of the secondary lattice and 
corresponding to the beating of the two lattices. The position of the secondary peaks is at Q = ±— (1 — 0). In the 
localised regime, the central peak of the momentum distribution is broadened till becomes of the order of the Brillouin 
zone (Fig. [8] (a) and (b))), while the secondary peaks are strongly suppressed. This is in agreement the experimental 
findings 15]. The broadening of the distribution in momentum space indicates that the states in real space are more 
and more localised. As interactions are turned on, we see a restoration of the peaks of the momentum distribution as 
a signature of the enhanced superfluid behaviour found by the analysis of the superfluid fraction. 

The transition from extended to localised phase has been characterised in the experiment [l5j by studying the 
visibility of the secondary peaks in the momentum distribution. In Fig. [5] we show the visibility of the momentum 
distribution, defined as V = ]n(0) — n(fci)]/[n(0) + n(kx)], both for the noninteracting and for the interacting case. 
The loss of visibility starts to become important at values of the height of the secondary lattice A/t which depend on 
the choice of (3 and of the interaction strength. In the case of an irrational (3 = (\/5 — l)/2 and in the noninteracting 
limit the visibility starts to decrease right at the transition point A/t = 2, while for other choices of (3 as in the 
experiments the visibility remains close to one to larger values of A/t, indicating a smoother crossover from extended 
to localised states. Similarly, in the interacting case the visibility starts to decrease from its value one at larger value 
of A/t, reflecting the shift of the transition point. 



10 



10 







c3 



10 



-5 




t s 



N=13 
N=21 
N=34 



10 







NIA-A 



v 



10' 



FIG. 7: (Colour online) Double logarithmic scaling plot of the superfluid fraction in the case U/t = 0.02. The parameters are 
A c = 3.2t, v = 1.25, a = 0.38. 



V. SUMMARY AND CONCLUDING REMARKS 



In summary, we have analysed the superfluid to Anderson-localisation transition for an interacting Bose gas sub- 
jected to a quasiperiodic potential, focussing on the case of a finite-size lattice. We have investigated the effects of 
changing the values of the quasiperiodicity parameter (3 and of the interactions. In the case of periodic boundary 
conditions, by studying the superfluid fraction we have found that the transition is sharper for irrational values of /3 
and smoother when several minima in the potential have almost the same energy. We have also found that the transi- 
tion is shifted to higher values of the secondary lattice height when the repulsive interactions are taken into account. 
Implementation of this idea could lead to an accurate method to measure of the interaction coupling strength, which 
is directly related to the s-wave scattering length . In the case of open boundary conditions, the inverse participation 
ratio is introduced as a good indicator of the transition, and its scaling analysis reproduces again the change of char- 
acter of the transition from irrational values of (3 to rational ones (with a large denominator). We have furthermore 
shown that also in the interacting case the change of the ground state of the system might be probed by analysing the 
momentum distribution of the gas, which displays side peaks in the superfluid phase but not in the deeply localised 
phase, and the visibility of the momentum distribution. In perspective, a detailed analysis of the shift of the transition 
at increasing the interaction strength, as well as the inclusion of a harmonically trapping potential [4||, seem to us 
valuable extensions of the present work. 
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FIG. 8: (Colour online) Momentum distribution (in units of the primary lattice spacing a) of a Bose gas in the quasiperiodic 
lattice, for various choices of the height of the secondary lattice A/t = 1 (solid lines), A/t = 2 (dashed lines) and A/t = 3 
(dot-dashed lines); of the interaction strength U/t = (left panels) and U/t — 0.02 (right panels) and of the periodicity 
parameter (3 — (\/E — l)/2 (top panels) and 1032/862 (bottom panels). For all the curves we have used a number of lattice 
sites N — 50 and a filling factor 1/2. 
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FIG. 9: Visibility of the momentum distribution (adimensional) of a Bose gas in the quasiperiodic lattice, as a function of the 
height of the secondary lattice A/t for three choices of the quasiperiodicity parameter (3 (solid line: (3 — (s/E — l)/2, dotted 
line j3 — 830/1076, dash-dotted line /3 = 1032/862); left panel, noninteracting case; right panel interacting case with U = 0.02t. 
For all the curves we have used a number of lattice sites JV = 50 and a filling factor 1/2. 
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